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Polariton polarization can be described in terms of a pseudospin which can be oriented along 
the X, y, or z axis, similarly to electron and hole spin. Unlike electrons and holes where time- 
reversal symmetry requires that the spin-orbit interaction be odd in the momentum, the analogue 
of the spin-orbit interaction for polaritons, the so-called TE-TM splitting, is even in the momentum. 
We calculate and compare spin transport of polariton, electron, and hole systems, in the diffusive 
regime of many scatterings. After dimensional rescaling diffusive systems with spatially uniform 
particle densities have identical dynamics, regardless of the particle type. Differences between the 
three particles appear in spatially non-uniform systems, with pumps at a specific localized point. 
We consider both oscillating pumps and transient (delta-function) pumps. In such systems each 
particle type produces distinctive spin patterns. The particles can be distinguished by their differing 
spatial multipole character, their response and resonances in a perpendicular magnetic held, and 
their relative magnitude which is largest for electrons and weakest for holes. These patterns are 
manifested both in response to unpolarized pumps which produce in-plane and perpendicular spin 
signals, and to polarized pumps where the spin precesses from in-plane to out-of-plane and vice 
versa. These results will be useful for designing systems with large spin polarization signals, for 
identifying the dominant spin-orbit interaction and measuring subdominant terms in experimental 
devices, and for measuring the scattering time and the spin-orbit coupling’s magnitude. 

PACS numbers: 71.36.+C, 72.25.Fe, 72.25.Dc 


I. INTRODUCTION 

Quantum particles are characterized not only by 
their position, but also by additional quantum num¬ 
bers. Among these quantities the electron’s spin, a prop¬ 
erty which has one of two values: +1/2, and —1/2, 
has tremendous promise for improvements in computing 
technology. If there is no magnetic field which breaks 
time reversal symmetry, then a 180 degree rotation re¬ 
versing an electron’s momentum must also reverse the 
evolution of its spin. Here we will explore diffusive spin 
transport of an alternate particle, which like electrons has 
a quantum number with two possible values, but which 
unlike electrons keeps the same dynamics when momen¬ 
tum is reversed. 

The particle of interest is the cavity polariton, the ele¬ 
mentary excitation of semiconductor microcavities in the 
limit of strong coupling between light and the quantum 
well’s excitons. Like electron spin, the polariton has a 
quantum number with two states: its polarization. Un¬ 
like electrons, reversing a polariton’s momentum has no 
effect at all on its polarization. To reverse the polar¬ 
ization’s dynamics one must rotate the polariton by 90 
degrees not 180 degrees, and thus interchange the x and 
y axis. (For a review of polarization properties of polari¬ 
tons see e.g. Ref. m) The consequences of this profound 
difference for polarization transmission under repeated 
scattering are at the heart of the current article. 

The electron’s single particle density matrix p is 
characterized by charge density JV and spin densities 


Sx, Sy, Sz and can be represented as a vector 

P = [V, 5,,, Sy, Sz] . (1) 

Like the electron, the polariton also has a density matrix 
composed of its particle density N, circular polarization 
density Sz, and linear polarization densities Sx, Sy cor¬ 
responding to xy and diagonal in-plane polarization. In 
the ballistic regime, i.e without scattering, the polaritonic 
analog of spin-orbit coupling causes the polarization’s di¬ 
rection to precess, converting linear Sx,Sy polarization 
to circular Sz polarization and vice versa. For polaritons 
the axis of precession varies continuously over 360 degrees 
when the polariton momentum rotates by 180 degrees, 
unlike electrons where the axis of precession is locked to 
the momentum. Because of this angular structure, an 
initial population of linearly polarized Sx or Sy polari¬ 
tons will evolve into into a four-lobed, or quadrupole, 
pattern of circular Sz polarization. In this celebrated 
effect, two lobes have positive circular polarization, two 
have negative polarization, and separating the four lobes 
is the Langbein cross where the initial linear polarization 
is unchanged. ^2H11 When a perpendicular magnetic field is 
added, the fo ur-q uadrant polarization pattern twists into 
spiral shapes.l^ 

We will compare the polariton’s dynamics under scat¬ 
tering to those of electrons, and also to holes realized 
in zinc blende semiconductors, where the Fermi level is 
doped into the J = 3/2 valence band rather than the 
conduction band. In each of these cases we will restrict 
our analysis to two-dimensional systems analogous to the 
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2-D electron gas. Like polaritons, electrons and holes 
can be written mathematically as two-state systems, and 
their dynamics are governed by two-state Hamiltonians. 
However, unlike polaritons, they possess real spin and 
therefore must reverse their spin dynamics when their 
momentum is rotated by 180 degre es. Moreover, in holes 
(without strain or anisotropjffSHil) rotating the momen¬ 
tum by 60 degrees also reverses the spin dynamics. This 
is a consequence of the crystal symmetry, in combina¬ 
tion with the four-fold degeneracy of the valence band.^^ 
Therefore an initial population of holes will evolve 
into a six-lobed, or sextupole, pattern of Sz spin density. 
Holes produce six lobes, polaritons produce four lobes, 
and electrons produce a simple dipole pattern. These 
qualitative differences in angular dependence correspond 
to precise differences in the two-state Hamiltonians gov¬ 
erning electrons, polaritons, and holes. 

In addition to the distinctive spin-orbit splitting terms 
of each type of particle, we will consider also the effect 
of a Zeeman splitting term caused by a external mag¬ 
netic field oriented perpendicularly to the sample. For 
polaritons this term may also be caused by variations in 
the polariton cavity such as in-plane strains inside the 
quantum well, or asymmetries in the direction of crys¬ 
tal growth.® This external and tunable parameter is a 
useful probe of (pseudo)spin dynamics, and in polariton 
systems is known to twist the polarization pattern into a 
spiral. 

Previous polariton experiments have generally utilized 
GaAs samples where the clean sample’s low disorder den¬ 
sity is further diminished by the polariton’s large spa¬ 
tial extent, allowing polaritons to move ballistically, i.e. 
without scattering, on millimeter length scales.^^^l*^ Bal¬ 
listic polariton transport has been studied both experi¬ 
mentally and theoretically in many papers.®® It is com¬ 
pletely controlled by (pseudo)spin precession: As the po¬ 
laritons move radially out from the polariton pump, their 
quadrupole pattern changes sign twice for each preces¬ 
sion length Iprcsn- The result is that concentric circles 
of alternating sign are superimposed on the quadrupole 
pattern.® Perfect spin polarization is maintained until the 
polaritons decay. 

In this article we focus on a different scenario, the dif¬ 
fusive regime, where the scattering length is so short that 
transport is governed by diffusion. We assume that the 
experimental length and time scales are long compared 
to the scattering scale. This is the typical case for experi¬ 
ments on 2-D electron and hole gases. For polaritons, the 
diffusive regime can be obtained by fabricating nanopil¬ 
lars embedded in the polariton sample, similarly to opti¬ 
cal systems where disorder has been induced by growing 
nanocolumns.*^^*^ In the diffusive case the particle’s av¬ 
erage displacement tq grows according to the diffusive 
law 

ro = V~Dt (2) 

where D is the diffusion constant, unlike the ballistic 
regime where the rms value of ro grows linearly with time. 


In diffusive samples the momentum becomes randomized 
through elastic scattering, subject to the constraint that 
the energy remains equal to the pump energy Ep. There¬ 
fore one calculates the average dynamics and transport - 
averaging over the random disorder, and averaging over 
the randomized momentum. Performing this averaging, 
and assuming a low-density, non-interacting system, we 
calculate the spin diffusion equations which govern the 
evolution of spin and particle densities. Previous works 
on polaritons in the diffusive regime used a classical ki¬ 
netic equation to calculate the e quilibrium of steady-state 
spatially uniform systems.l^SEIl Here we calculate the dy¬ 
namics of non-uniform systems using standard perturba¬ 
tive techniques based in quantum mechanics which have 
been applied extensively to electron diffusion. 

We will show that in the diffusive regime spatially uni¬ 
form distributions of electrons, polaritons, and holes ex¬ 
hibit identical spin dynamics. The essential reason is 
that the three particles differ only in their spin-orbit in¬ 
teraction’s dependence on momentum, and momentum 
is randomized by scattering in the diffusive regime. Dif¬ 
ferences between the spin diffusion equations of the three 
particles can be resolved only by breaking translational 
invariance, so we will focus entirely on non-uniform sys¬ 
tems. 

Previous articles deriving spin diffusion equations have 
often applied the equations to quasi-I-D g eomet ries in 
steady-state, i.e. without temporal dynamics 1221231 j 
case the diffusion equations simplify and can be solved 
analytically. In this present article we focus instead on 
the richer and physically relevant case of spatially local¬ 
ized pumps introducing particles into the device at one 
particular point. Such pumps produce a rich set of inter¬ 
esting spin behaviors. We will show that in the presence 
of a magnetic field, unpolarized pumps produce through 
scattering both in-plane and out-of-plane spin densities. 
Polarized pumps do the same, and in addition exhibit 
spin precession between in-plane and out-of-plane spins. 
These signals display spatial multipole patterns with the 
same angular structure as the spin-orbit coupling - a 
quadrupole structure in the case of polaritons. If there 
is a perpendicular magnetic field, the multipoles twist 
into spirals as time progresses. However, unlike ballistic 
transport, the spin pattern does not manifest concentric 
circles around the pump, since transport occurs accord¬ 
ing to a random walk rather than radial motion. 

The structure of this paper is as follows. In section [lT| 
we explain in detail the electron, polariton, and hole sys¬ 
tems which can be realized experimentally, and realistic 
Hamiltonians for modeling these systems. We also de¬ 
scribe how to quantitatively compare diffusive transport 
of these disparate particles, by rescaling the dimensions 
and comparing dimensionless ratios. In section in we 
mathematically prescribe the observables which are cal¬ 
culated in our paper - the density matrix which encodes 
spin and particle densities, and the response function for 
a spatially localized pump. Next section |IV] presents our 
analytical results: the diffusion equations which control 
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the evolution of spin densities. These diffusion equations 
are obtained using a standard technique which has been 
widely applied to electron diffusion problemsExam¬ 
ination of the diffusion equations reveals the multipole 
pattern which will allow experimentalists to quickly de¬ 
termine the dominant spin-orbit coupling and carefully 
measure sub-dominant couplings. Changing the spin- 
orbit splitting or the scattering time will cause the spa¬ 
tial multipole pattern to rotate. We also show that in a 
perpendicular magnetic field polaritons and holes exhibit 
a resonance that can be used to measure the spin-orbit 
strength. The resonance is characterized by extinction of 
the polarization produced in response to an unpolarized 
pump. It occurs at numerically different magnetic field 
strengths for the two particles, and is sensitive to the spin 
splitting. 

The behavior of a real experimental device can not be 


determined unless the diffusion equations of section IV 


are combined with initial conditions and solved. Unlike 
previous works which analytically calculated steady-state 
quasi-1-D devices, in section|V]we solve 2-D spin diffusion 
from a localized pump, using numerical techniques. The 
problem/solution space is very rich: oscillating pumps vs. 
transient pumps, four types of spin input and response 
{Sx, Sy, Sz, and density N), various angular and radial 
profiles of the response, time dependence at both short 
and long times, and an intricate dependence on the rel¬ 
ative strengths of the kinetic energy, spin-orbit splitting, 
Zeeman field, and scattering time. Surveying this com¬ 
plex solution space was one of the chief challenges of the 
present work. Section |V] presents the most interesting 
features. In many cases our results quantitatively con¬ 
firm and illustrate the qualitative trends discussed in the 
previous section W However Sect ion |V] also gives qualita¬ 
tively new information: we show that the time evolution 
breaks naturally into two regimes, before and after loss of 
memory of the pump polarization. Spin observables can 
reach quite large values in the first regime, and can per¬ 
sist well beyond the memory loss. We also examine the 
evolution into spiral structures which will allow measure¬ 
ment of the scattering time, and determine the param¬ 
eters which maximize the spin polarization. Lastly Sec¬ 
tion |VI| synthesizes the major differences between elec¬ 
trons, polaritons, and holes. Several appendices present 
mathematical details of our calculations. 


II. THE MICROSCOPIC HAMILTONIANS FOR 
ELECTRONS, POLARITONS, AND HOLES 


We now introduce a realistic Hamiltonian which mod¬ 
els polaritons in a GaAs quantum well, as well as elec¬ 
tronic 2DEGs and hole-doped systems. The spin-orbit 
coupling varies substantially in these three systems. It 
is the linear Rashba coupling for electrons, a quadratic 
coupling for polaritons, and a cubic coupling for holes. 
We write the Hamiltonian for all three cases in a unified 
notation with TV = I for electrons, TV = 2 for polaritons. 


and TV = 3 for holes: 
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where tf = [o'x,o'y,az] are the Pauli sigma matrices, k is 
the wave vector, and kp is the characteristic wave-vector 
determined by the energy E. For convenience we have ex¬ 
pressed the relative strength of the Zeeman term and the 
spin-orbit coupling with the angle Op = arctan(Afc^ /bz), 
a sine S = sin 0 b, and a cosine C = cos 0 b- ‘2-Eprcsn = 

2 ^(Afcp )2 -I- 62 jg splitting between energy levels. 

This is a non-interacting Hamiltonian, which works 
well for polaritons in the limit of low polariton density, 
where polariton interactions are weak and can be ne¬ 
glected. In this limit the polariton energy Ep - the in¬ 
plane kinetic energy in the lower polariton branch - is 
the same as the energy of the pump photons used to ex¬ 
cite the polaritons. Moreover the polaritons’ character¬ 
istic wavelength is the magnitude of k on the contour 
in momentum space defined by E{k) = Ep. (This con¬ 
tour is the elastic scattering circle for polaritons, or the 
Fermi surface for electrons and holes.) Note that in the 
regime of strong pumps nonlinear effects m ay qu alita- 
tively change the texture of the spin patternPIESl xhe 
non-interacting picture is also appropriate for electrons 
in 2DFGs and for hole-doped semiconductors, with the 
caveat that here Ep is the Fermi energy Ep, and is con¬ 
trolled by doping, gating, etc. 

Our TV = 2 polariton Hamiltonian is accurate when 
the polariton energy Ep is a fraction of Rabi splitting H-b, 
which is the energy scale of the exciton - photon coupling 
that produces the polaritons. In GaAs the Rabi splitting 
is near the value = 6meV seen in Refs. mHHEnsni 
so our Hamiltonian is valid in the range Ep = [0, 2] meV. 
In this article we will fix Ep at 2.0 meU; decreasing this 
value only rescales the scattering length and time, as well 
as the polariton wavelength Ip determined by the pump 
energy. We also choose the mass m = 5 x 10“^me, as 
reported in References mHT51 This value is typical for 
polaritons in GaAs quantum wells, and is determined by 
the band gap in GaAs. 

The polarization evolution (or spin evolution for elec¬ 
trons and holes) is determined by the second and third 
terms in H. The second term is a Zeeman splitting de¬ 
pending on the polarization (or spin) quantum number. 
Here we assume that this term is oriented along the z 
axis, perpendicularly to the sample, so it produces only 
in-plane precession. 
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The third term gives the momentum-dependent spin- 
orbit coupling, or the TE-TM splitting in the case of the 
cavity polariton. In polariton systems symmetry argu¬ 
ments require that this term be even in the momentum, 
which is a very distinctive signature of polaritons. In 
contrast, in electron and hole systems time reversal sym¬ 
metry requires that this term be odd in the momentum, 
as long as there is no magnetic held or magnetic impuri¬ 
ties. The = I linear term we use here for electrons is 
just the Rashba spin-orbit interaction. 

The principal source of the quadratic polariton TE- 
TM splitting is the conhned photons’ sensitivity to the 
mirrors of the cavity and to the angle of the incident 
light.® Studies of the TE-TM term and of the Zeeman 
term (with a modest 5 Tesla held) have shown that both 
terms can reach or exceed 0.2 meE, so that they may be 
comparable to the total polar iton energy, and spin-orbit 
effects can be very large 


A. How to Compare Electrons, Polaritons, and 
Holes 

For the purpose of comparison between uncharged po¬ 
laritons and charged electrons and holes, we neglect the 
electrical potential $ and magnetic gauge potential A 
acting on charges. The gauge potential A can be ne¬ 
glected if the cyclotron radius is larger than both the 
scattering length and the dephasing length, which is eas¬ 
ily achieved if the sample is not very clean and the tem¬ 
perature is high enough to extinguish weak localization 
effects. Turning to the electrical potential, a gradient in 
this potential will cause a net drift of the charges, which 
will act in addition to the diffusive dynamics considered 
here. This drift can be be minimized by decreasing the 
potential gradient. 

In the preceding discussion of the Hamiltonian we have 
taken care to use realistic values for the polariton mass to, 
pump energy Bp, and spin-orbit coupling strength. How¬ 
ever the final results of our calculations are determined 
only by length and time scales and dimensionless ratios 
formed from combinations of these parameters. We can 
compare polariton results directly to electron and hole 
results by rescaling the dimensions to match each other. 
Here are the physically relevant scales: 

The Wavelength Ip and Velocity Vp. The mass to 
and the pump energy Ep (or Fermi energy for electrons 
and holes) determine the wavelength 

Ip = h/y/2mEp, (4) 

which is the fundamental length scale governing trans¬ 
port and dynamics. They also determine the polariton’s 
velocity in ballistic flight, 

Vp = yj2Ep/m, (5) 

which conveniently allows conversion of times and ener¬ 
gies to lengths, and vice versa. At the polariton energy 


Ep = 2.0 meV the polariton wavelength and velocity are 
respectively Ip = 4^to and Vp = Afim/ps. 

The Precession Length Iprcsn- Taken together, the 
Zeeman and spin-orbit terms determine the spin (or po¬ 
larization) precession length 

iprcsn — ^'^p! ^prcsnj ( 6 ) 


where Eprcsn = \J+ b'^. We assume that the 
spin-orbit term is small compared to the kinetic energy 
/2m, or equivalently, that the precession length Iprcsn 
is large compared to Ip = A pm. 

The Scattering Length I scales with 


I oc 


{HveY 
(V2) y2 


( 7 ) 


where (V^) the second moment of the disorder potential 
and V is the characteristic length scale of the impurities. 
In the diffusive regime only I itself determines the sys¬ 
tem’s evolution, not n or (V^). We do require that the 
disorder be weak according to the loffe-Regel criterion; 
i.e. that x = Ip/lir ^ 1. Otherwise the diffusion pic¬ 
ture breaks down and Anderson-localized states may be 
observed. To ensure diffusive transport in our polariton 
system with Ip = A pm, we set I = 16 pm. The scattering 
time is T = l/vp = Aps. 

The Polariton Decay Length Idecay Unlike elec¬ 
tron and hole systems where charge is conserved, polari¬ 
ton number generically decreases with time as photons 
escape through the Bragg mirrors of the quantum cavity. 
This decay is determined by the polariton lifetime r^ecay, 
which also determines the maximum distance Idecay that 
polaritons travel before decaying. Under diffusive motion 
this distance is 


^decay — \/ decay 


( 8 ) 


where D = lvp/2 is the 2-D diffusion constant. Diffusion 
can be observed only if the decay length Idecay substan¬ 
tially exceeds the scattering length I = 16 pm. Therefore 
we assume a relatively large value Idecay = 80 pm, cor¬ 
responding to Tdecay = 200ps = 50t. The (pseudo)spin 
dynamics of greatest interest occur at shorter time scales 
of order lOr, so our numerical results on spin are largely 
independent of the decay time. For electrons and holes 
the decay time is infinite. 

These numerical values of the polariton decay length 
and time are experimentally accessible. Polariton cavities 
far larger than Idecay = 80 pm are routinely fabricated 
in present experiments. Improved mirrors have recently 
increased the lifetime from a few ps to 180 ps, so 200ps 
is also quite feasible 

Given these length scales, electrons, holes, and polari¬ 
tons can be compared directly to each other after posi¬ 
tions are normalized by the scattering length I and times 
are normalized by t. These are the natural length and 
time scales of a diffusive system. After performing this 
normalization, only four dimensionless quantities deter¬ 
mine the (pseudo)spin dynamics: 









5 


• The dimensionless disorder strength x = = 

h/EpT <c 1. 

• Ob, which controls the relative strength of the Zee- 
man and spin-orbit terms. 

• The dimensionless diffusion length, vq/I = ■\Jtl2T, 
which tracks the diffusive spreading of the initial 
particles. 

• The dimensionless energy splitting between the 
spin up and down states ^ = l/lprcsn = EprcsuT/h 
which determines the dominant source of spin 
relaxation.SilH^ 

This last parameter is key to the spin dynamics. C ^ 1/2 
puts the system in the Elliott-Yafet (EY) regime, where 
the (pseudo)spin processes many times before scatter¬ 
ing and momentum is tightly coupled to spin. In the 
EY regime the spin is randomized at every scattering 
event, so its relaxation time is close to the scattering 
time T. In contrast, C 1/2 puts the system in the 
D’yakonov-Perel’ (DP) regime where the (pseudo)spin 
precession length Iprcsn is long compared to /, and the 
energy splitting is small compared to the scattering en¬ 
ergy. In this case precession lasts much longer than the 
scattering time, and the spin relaxation length scales with 

iprcsn • 

As always, we assume that Ip/Iprcsn = ’’"Cx is small, so 
that the spin splitting is small compared to the kinetic 
energy p^/2m. 

III. THE DENSITY MATRIX AND THE 
RESPONSE FUNCTION 

We now begin developing the equations describing 
(pseudo)spin dynamics in the diffusive regime, where the 
particle scatters randomly many times, and where we 
must average over disorder in order to calculate transport 
properties. Disorder-averaged diffusive transport is qual¬ 
itatively different from transport in the ballistic regime, 
where all transport information is contained in the wave- 
function !'(/)). In the diffusive regime scattering random¬ 
izes |'0)’s phase and therefore its disorder average (|V'))e 
averages to zero, losing all transport information. Study 
of the disorder-averaged wave-function alone is unable to 
describe transport, or even probability conservation, in 
the diffusive regime. 

A correct description of diffusive transport begins with 
the single-particle density matrix which for the case of a 
pure state reads: 

pit) = l^(i)) ® (V'(^)l (9) 

pit,x,x) = {xli/it)) {-tp{t)\x) 

The density matrix encodes complete information about 
the probability density I{x) in its local x = x values: 

( 10 ) 


Here the trace is performed over p’s spin indices. I{x) is 
proportional to the charge density in electron and hole 
systems, and to the particle density in polariton systems. 
Similarly, the (pseudo)spin densities are also encoded in 
p’s local values: 

Sxit,x) = Tr{p{t,x,x)ax) 

Sy{t,x) = Tr{p{t,x,x) ay) 

Szit,x) = Tr{p{t,x,x)a^) (11) 

These local quantities survive the disorder average be¬ 
cause of the cancellation of phases between \'tp) and 
(■01. Nonlocal x ^ x values of the density matrix de¬ 
cay exponentially when a: — i? is longer than the scatter¬ 
ing length if we are outside the polariton condensation 
regime^SI. We confine our attention to the local quantities 
I{t,x), Sx/t^x), Sy(t,x), Sz{t,x), which are the proper 
starting point for studying diffusive transport. We write 
these as a vector: 

(p(t,£))'^ = {lit,x), Sxit.x), Sy{t,x), Sz{t,x)} {12) 

It is also convenient to study the spatial Fourier trans¬ 
form of the density matrix, because random disorder af¬ 
ter averaging does not break the (average) translational 
invariance: 

ipit,^)'^ = [lit,^, Sxit,^, Sy{t,^, Szit,^}il‘i) 

When studying transport with an oscillating pump we 
will also study the temporal Fourier transform 

{p{u},x))’^ = j du} exTp{iujt) {p{t,x))’^ (14) 

It is often useful to compare the spin densities Si to 
the probability density I. We define the (pseudo)spin 
degrees 

Sx = Sx/I, Sy = Sy/I, Sz = Sz/I (15) 

These degrees can never be larger than 1 or smaller than 
— 1. At one extreme, the density matrix may be en¬ 
tirely unpolarized, in which case it is proportional to 
p = {1,0,0,0} and the spin degrees are all zero. At 
the other extreme it may be completely polarized, with 
+ Sy + = 1- 

A. Response Functions 

Our goal is to obtain and study a response function 
$(f) which evolves the density matrix forward in time: 

rXi) = = (16) 

$(t, if, f) is a 4 X 4 matrix, because it computes the charge 
and spin densities at time t as function of the charge and 
spin densities at t = 0. The response function is also 


Iit,x) = Tr{p{t,x,x)) 
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a linear operator, which is a direct consequence of our 
working in the noninteracting regime. Therefore we will 
study only the response $(<, if) to a t = 0 delta function 
input po (x) at the origin. The response to any input 
may be reconstructed from ^(t,x) by convolving the re¬ 
sponse function over the t = 0 density matrix according 
to equation [THl 

We will also calculate the temporal Fourier transform 
<i>(a;,if) which describes the density matrix produced in 
response to a pump po oscillating with period T = 27r/a;: 

<i)(a;,if) = J dt exp{ujjt) $(t,if) 

p{u},x) = / dx ^{u},x — x) po{uj,x) (17) 


\ G^(k^-ql2) 

1 

1 G^(k 2 - q/2) 

1 

1 

A 

^ k-2 — 

G^(+g/2) 1 G^{kr+q/2) 

.1 ^_ 

1 G^{k 2 + q/2) 

1 — 


FIG. 1: Diagram illustrating the scattering that produces 
diffusion. Here two scattering events are shown. and 
describe time evolution of the polariton ip and its complex 
conjugate ip^. Each scattering event causes correlations be¬ 
tween G'^ and G^ and is shown as a dashed line connecting 
the two. 


This frequency response function contains both informa¬ 
tion about the response amplitude, and also information 
about the response’s phase relative to the pump. There¬ 
fore it is in general complex, with both real and imagi¬ 
nary components. In this paper we will study only the 
magnitude of the response, not its phase information. 

Because the system is translationally invariant after 
averaging over disorder, it is convenient to Fourier trans¬ 
form both the time response function $(t,af) and fre¬ 
quency response function <i>(a;,af) to momentum space: 

p{t,x) = J -^0^ exp(ig • f) = 0,^ (18) 

We determine the response functions using standard 
methods from the diagrammatic technique for disordered 
system^. We neglect effects from quantum interference 
and adopt a purely classical picture of particle diffusion 
through the disordered system. We model scattering 
with a ’’white noise” disorder potential which does not 
alter the (pseudo) spin 


where the density of states v{Ep) at pump energy Ep is 

^^Ep) = E/ (20) 

and E[s, k) are the two eigenvalues of the Hamiltonian 
at wave-vector k. We assume that this equilibrium value 
of p is unpolarized, which just means that Ep > 2Eprcsn- 
With this white noise disorder every scattering event 
causes the particle to lose its memory of its previous mo¬ 
mentum, which becomes evenly distributed on the elastic 
scattering circle. Under these assumptions the time re¬ 
sponse function can be written as an exponential: 

d>(t,^ = exp{-t > 0 (21) 

T 

D~^{q) = 1 — lij is the inverse of an operator called the 
diffuson. It contains information about a single scatter¬ 
ing event which is represented by the scattering operator 


u{f), {u{E)u{^) = {2TTVT/h) ^i5(r —-^(19) 


lij. The response function is obtained by exponentiat¬ 
ing describing the results of repeated scatterings. 

As a consequence of this exponential form, the frequency 
response function is 




1 

lU) + T~^D~^{q) 


( 22 ) 


Equivalently to equations and we can also write 
the spin diffusion equation 


= D ^{q)p{t,q). 


(23) 


This reveals that D~^ is the diffusion operator control¬ 
ling evolution of the density matrix. Before solving the 
diffusion equation one must first supplement it with the 
starting state p{t = 0). 

In the next section we will derive the diffusion operator 
D~^{q) analytically. Calculation of the response func¬ 
tions requires numerical calculations - either exponenti¬ 
ation (Equation |21[ ) or inversion (Equation \22\ - and we 
will perform this numerical work in Section]V[ 


IV. ANALYTICAL CALCULATION OF 
DIFFUSIVE TRANSPORT 

In order to calculate the diffusion operator one 

must first calculate the scattering operator lij. In each 
scattering event the wave-function \il}) and its conjugate 
('01 move together, scattering in unison. Two scattering 
events are pictured in Figure The scattering operator 
lij is given by the integral 

= 0^ J dkTr{G^ik-q/2,Ep)a, 

X G^{k + q/2,Ep)aj) (24) 

Here G^ and G^ are the disorder-averaged single-particle 
Green’s functions, and their spectral representation is 
given in equation |C2| in Appendix [Cj The trace is taken 
over the spin indices of G^, G^, ai, and aj, which are all 
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2x2 matrices in (pseudo)spin space, q is the diffuson 
wave-vector. Further details of the integration of equa¬ 
tion are given in equations |C3| and |C4| of Appendix 

o 

Equation!^ for the scattering operator, and = 

1 — lij for the diffuson, are a well established and widely 
adopted formalism for determining spin and charge 
diffusion.l22H2n! This formalism has a strong physical mo¬ 
tivation expressed in the perturbative diagram in Fig¬ 
ure [2 can be derived from the Keldysh equations for 
nonequilibrium conduction, and is equivalent to Kubo 
linear response theory.l^ In the present work we use this 
formalism to study electrons, polaritons, and holes, in a 
magnetic held, to all orders of the spin splitting param¬ 
eter (. We alsos have checked that our results coincide 
with those of Reference [22l which is also based on the 
Keldysh equations, but uses a different strategy for solv¬ 
ing them iteratively. 

As seen from the Fourier transform in equation |18[ 
there is a direct mapping from the momentum q to posi¬ 
tion space. Precisely, every power of tqx corresponds to 
a d/dx spatial derivative, every power of iqy corresponds 
to a d/dy, and %q corresponds to V. Because we are in¬ 
terested only in length scales that are long compared to 
the scattering length I, we perform a Taylor series ex¬ 
pansion of the scattering operator lij in powers of Iq (or 
equivalently in powers of ZV) and terminate the expan¬ 
sion at second order. This approximation is justihed for 
calculating observables at scales larger than 1. It guaran¬ 
tees that the inverse diffuson contains only terms with at 
most two spatial derivatives, which is appropriate fo r de- 
scribing particle diffusion under random scattering.l22l2ll 
For the spin-charge coupling in holes we make an excep¬ 
tion to this procedure and keep cubic (Iq)^ terms, since 
at quadratic order thi s coupli ng is zero, as is well known 
from previous studies 

As mentioned before, we assume that disorder is weak, 
far from the localized phase, so that the dimensionless 
disorder strength X = Z^/ttZ = h/EpT is small. Therefore 
we expand lij in a Taylor series in powers of x trun¬ 
cate at first order. We neglect all powers of x that do 
not occur in combination with Z, r, or C.. 

Our spin-orbit term has the special property that the 

spin splitting 2Eprcsn = 2^/(AZc^)^ -I- bl is rotationally 
invariant. As a result we are able to analytically calculate 
the scattering operator lij to all orders in the dimension¬ 
less energy splitting C = l/lprcsn- 


A. The diffuson in the EY and DP regimes 

After Taylor expanding equation and integrating 
over the elastic scattering circle, we find that the in¬ 
verse diffuson D~^ decomposes into four parts: a mo¬ 
mentum independent part A which gives the spin (or po¬ 
larization) relaxation rates and the precession term, a 
diffusion term A which includes the coupling between 


charge/number density and circular Sz polarization, a 
momentum-dependent part kq; which gives the coupling 
between density and linear Sx,Sy polarization, and a 
momentum-dependent part kic which gives the coupling 
between linear Sx , Sy polarization and circular Sz polar¬ 
ization. Our decomposition reads: 


D — (A— A-f Koi -|- Kic + T I'^'decay) (25) 


We present the four terms both in the Elliott-Yafet limit 
where C ^ 1/2, and also in the D’yakonov-Perel’ limit 
where <C 1/2. In Appendixj^we present the more com¬ 
plicated expressions which interpolate between these lim¬ 
its. Appendix [C| shows gives the details of our derivation 
of kqi and also of the coupling between particle density 
and Sz polarization. 

Spatially uniform spin distributions are governed by 
the mass and precession matrix. Very remarkably, this 
term is the same for electron, polariton, and hole-doped 
systems, regardless oi N = 1,2,3. It is: 






■0 0 0 0 - 

0 2(2(1-fC2) 2CC 0 

0 -2CC' 2(2(1-bC^) 0 

.0 0 0 4 ( 25 ' 2 . 


■0 0 0 0 ■ 

0 (l-bC'2)/2 0 0 

0 0 (l-bC'2)/2 0 

0 0 0 52 


(26) 


where the C and S parameters are defined in Eq. 

The on-diagonal terms in A give the (pseudo)spin re¬ 
laxation rates, which are the inverses of the spin life¬ 
times. The upper left entry is the inverse lifetime of 
the charge/number density, and is zero because elastic 
scattering does not change the polariton number or the 
electron or hole charge. In the EY ( 3> 1/2 regime the 
lifetimes are proportional to the scattering time r = 4ps, 
while in the DP ( <C 1/2 regime they scale quadratically 
with (“^; Txy.Tz oc r/(2 oc . 

The off-diagonal terms of A (zero in the EY regime) 
describe the usual precession of the in-plane polarization 
around the magnetic field, which is oriented perpendic¬ 
ularly to the sample, along the z axis. They render the 
diffuson non-Hermitian, so that its eigenvalues are com¬ 
plex, as required for describing precession. In the DP 
regime the polarization processes many - 0((“^) - times 
before it decays. In contrast, in the EY regime polariza¬ 
tion precession is not visible at diffusive scales, because 
each scattering randomizes the polarization. 

We turn to A, the part of the diffuson which is pro¬ 
portional to q2 and therefore describes diffusion. Since 
it is zero at zero momentum g = 0, it has no effect on 











spatially uniform distributions. 


- 


doz,N 



r 1 

0 

0 

doz 

{qir 

0 

1 

-6CC 

0 

2 

0 

6CC 

1 

0 


doz 

0 

0 

1 J 


r 1 

0 

0 

doz' 

(qir 

0 

5-2/2 

0 

0 

2 

0 

0 

52/2 

0 


doz 

0 

0 

C2j 


-XCC + xCCS^ N 

XC/{1 + N (27) 


The most interesting aspect of this diffusion operator 
is (Iqz, the coupling which controls the generation of Sz 
circular polarization from polariton density (or electron 
or hole charge). This coupling implies that from an initial 
population of unpolarized polaritons, diffusive scattering 
in a magnetic field will produce an Sz signal perpendicu¬ 
lar to the plane, with a conversion efficiency controlled by 
xC = Ip/T^lprcsn- The resulting Sz signal is rotationally 
symmetric in real space. The doz coupling describes dif¬ 
fusive generation of Sz polarization inside the structure 
of the diffuson D~^ itself. It is not the same as the exter¬ 
nal Zeeman potential introduced in Ref. [2H which does 
not couple Sz to density, and is not part of the diffuson. 
The doz coupling requires inclusion of the Zeeman split¬ 
ting while computing D~^] it is zero when — 0, and 
therefore is omitte d in s ome previous works on diffusion 
with a Rasha term .1221^ 

doz’s magnitude xCid^Y the same for electrons, po¬ 
laritons, and holes, but its dependence on 9b, the bal¬ 
ance between Zeeman and spin-orbit terms, differs for 
the three particles. We focus on the EY regime where 
the last term in d^z. can be neglected. In an electron sys¬ 
tem with linear spin-orbit coupling d^z does not depend 
on the spin-orbit strength A at all. Interestingly, in po¬ 
lariton systems if the magnetic held strength \ bz \ = \^kp\ 
is tuned to exact resonance with the spin-orbit strength, 
then doz goes to zero and changes sign. In hole-doped 
systems with cubic coupling the same resonance is still 
present, but is shifted to \bz\ = ^/2\Akp\. This reso¬ 
nant behavior is interesting both because it is different 
for electrons, polaritons, and holes, and also because it 
may be useful for measuring the spin-orbit strength. 

The diffusive component A contains an additional term 
A^ which is proportional to . Interestingly, it 

breaks rotational symmetry, mirroring the spin-orbit cou¬ 
pling. In electron systems with a linear Rashba spin-orbit 
coupling, A^ is: 


Ai 


ro 0 0 O' 

_ (^ 0 cos(26»,) sin(26»,) 0 

® 2 0 sin(2dg) — cos{26q) 0 ’ 

0 0 0 0 . 

2C2S'2(3 + 6C2-k8C^) 


(28) 


In the DP C <C 1/2 regime this term is proportional to 
which is a small contribution to A^^. However in the EY 
C 3> 1/2 regime where the spin-orbit coupling is strong 
A^ is proportional to S^/A, and induces a quadrupole an¬ 
gular dependence in the diffusion. A similar anisotropic 
diffusion term has been found for spin conduction on the 
surface of a 3-D TI, which be cause of its strong spin-orbit 
coupling is in the EY regime.^2MiJ More generally, in po¬ 
laritons A^ scales with and in holes it scales with 
(gZ)®. Because in polaritons and holes A^ is parametri¬ 
cally smaller than the (ql)^ part of A, we neglect it in 
these systems, and retain it only for electrons. 

Unlike Ref. |25Js treatment of the Rashba coupling 
which leaves the diffusion term constant at its ^ = 0 
value where the spin splitting is very small, our result for 
A keeps all orders of This allows spin or pseudospin to 
diffuse at a rate different from that of charge or number 
density, a possibility discussed in references [52HS41 While 

/ 7x2 

the diffusion term controlling the number/charge 
density is insensitive to (/, the other terms in A are re¬ 
duced by the spin-orbit coupling. Eor instance, the Sz 
diffusion term (expanded to second order in C) is reduced 
to ^^^(1 — 24(/^). This implies that spatial inhomo¬ 
geneities of number/charge density are smoothed more 
quickly than spatial inhomogeneities of (pseudo)spin. 

We find also a coupling which generates Sx, Sy linear 
polarization from polariton number density (or electron 
or hole charge density): 


Koi 


Xn 


cos N9q sin N9q 

0 0 

0 0 

0 0 

xC^ X {(I + C^)/A, CV2, (3/I6)(2C2 - 
xC/il + A^)^ S{l + C^) 
x{-1/4,-1/2,-(9/I6)(I-4CV3)} 


-{iql) XN 


0 

cos N9q 
sin N9q 
0 


(29) 


5^)} 


(30) 


This coupling implies that an initial population of unpo¬ 
larized particles, in the presence of a spin-orbit (TE-TM) 
splitting, will after diffusion produce a linearly polarized 
Sx, Sy signal with a conversion efficiency controlled by 
xC = Ip/'^lprcsn- This process does not occur without a 
spin-orbit coupling. 

Production of linear Sx,Sy polarization from an un¬ 
polarized source depends very sensitively on whether the 
system hosts electrons, polaritons, or holes. First of all, 
the polarization has a special angular pattern: a dipole 
pattern in an electron system, a quadrupole pattern in a 
polariton system, and a sextupole pattern in a hole-doped 
system. Secondly, the signal’s magnitude is strongly sen¬ 
sitive to the spin-orbit interaction, scaling with ql to the 
N — th power. At the length scales of interest, much 
larger than the scattering length I, this implies that the 
linear polarization is strongest for electrons and weak¬ 
est for holes. Thirdly, in the DP regime spin production 
from an unpolarized source is much smaller for electrons 
than for polaritons or holes. For electrons the two terms 
in equation nearly cancel, giving a signal that scales 


e 


(I + 4(2)3 
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as In contrast, for polaritons and holes there is no 
such cancellation and the signal is proportional to xC^ 
i.e. larger by two factors of the inverse dimensionless 
spin splitting which is large in the DP regime. 

Lastly, the linear polarization’s dependence on the Zee- 
man term is different for each of the three particles. The 
difference between particles is particularly strong in the 
EY regime where the second term in xn can be neglected, 
and relatively weak in the DP regime. We focus on the 
EY regime in the following discussion. In electrons the 
Zeeman term’s influence is weak, changing the polariza¬ 
tion by at most a factor of 2. In contrast, in polaritons 
removing the Zeeman term (setting (7 = 0) extinguishes 
the linear polarization. This means that polaritons, un¬ 
like electrons and holes, require a magnetic field to pro¬ 


duce a linear spin polarization signal from an unpolarized 
pump. In holes the requirement is again absent, but one 
can tune the magnetic field to zero the linear polarization 
at the resonant condition \bz\ = |A|/-\/2. As one moves 
from the EY regime to the DP regime the resonance shifts 
to smaller \bz \ and eventually disappears. This resonance 
may be useful for measuring the spin-orbit strength in 
hole-doped systems. 

Finally we turn to the coupling between linear S^, Sy 
polarization and circular Sz polarization, kic, which de¬ 
scribes out-of-plane (pseudo)spin precession. Unlike the 
couplings doz, kq/ to the number or charge, the form of kic 
depends strongly on the dimensionless energy splitting 
and is independent of the disorder strength x- It changes 
if one goes from the DP regime to the EY regime: 


Ulc 


DP : 
EY : 


■0 0 0 

/ i\N 0 0 0 

0 0 0 

.0 /at cos N9q-\-gifSinN9q Jn sin N9q — cos N9q 

fN=0, gAr = C^x{2,-3,l} 

fN = -CS{-2)-^, gN = 0 


0 

/at cos N9q — gN sin N9q 
/at sin N9q + gN cos N9q 
0 


(31) 


This coupling implies that a spin-orbit coupling will cause 
out-of-plane precession, which converts an initially lin¬ 
early polarized Sx, Sy signal to circular Sz (pseudo) spin, 
and vice versa. In the DP regime only a spin-orbit cou¬ 
pling is required, while in the EY regime a magnetic field 
is also required. In the DP regime this is bona fide spin 
precession, since the g^ contribution is anti-Hermitian, 
similarly to the precession term coupling Sx and Sy. In 
this regime oscillations between Sx, Sy and Sz are pos¬ 
sible. In the EY regime oscillations are not possible be¬ 
cause the /at term is dominant and is Hermitian; we find 
only a non-oscillatory coupling. 

Similarly to kqi which produces linear polarization 
from an unpolarized pump, this precession between lin¬ 
ear and circular polarization is strongest for electrons and 
weakest for holes, and the resulting polarization shows a 
dipole, quadrupole, or sextupole pattern for electrons, 
polaritons, and holes respectively. An interesting differ¬ 
ence is that the angular pattern rotates by tt/ {2N) during 
the transition from small C to large This angular pat¬ 
tern is expressed in the coupling’s dependence on sin N9q 
and cos N9q. These sines and cosines are interchanged 
when C is increased from the DP limit to the EY limit, 
causing the rotation. 


B. Instability and Regnlarization 

Our calculation of the response function involved a 
Taylor series expansion in powers of ql. This is a long- 


wavelength approximation, restricted to length scales 
larger than the scattering length. Outside of this ap¬ 
proximation’s range of validity, i.e. at momenta ql > 1 
large enough to probe scales smaller than the scattering 
length, our analytical formulas predict that the Sx,Sy 
linear polarization grows exponentially with time instead 
of decaying, which is an unphysical result. In polariton 
systems this instability occurs only when C ex 1/2 is near 
the transition between the EY and DP regimes, i.e. when 
the spin precession length matches the scattering length. 

We have regularized this instability by make an ad-hoc 
modification of the diffuson, i.e. a cutoff, when I < {qlY' . 
This cutoff is technically sound and justified because we 
already removed all short-distan ce phy sics when we trun¬ 
cated higher spatial derivatives We have recalcu¬ 

lated our numerical results with three different cutoffs, 
and we compared the results. When the computed re¬ 
sults differ, this signifies that they are sensitive to ballis¬ 
tic physics, i.e. length scales smaller than the scattering 
length 1. Throughout our numerical discussion we report 
which results are stable, independent of the cutoff, and 
which results do vary with the cutoff because of a sensi¬ 
tivity to ballistic physics. Appendix [B] describes each of 
the three cutoffs in more detail. 

The numerical results presented in our graphs are cal¬ 
culated with a smooth cutoff that turns on smoothly over 
the interval | < {ql)^ < 2. This cutoff nonetheless pro¬ 
duces a spatial oscillation, or ringing, in some of our re¬ 
sults. The ringing has a length scale set by the cutoff 
scale I and affects all components of the density matrix. 
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C. Summary of Analytical Results 

Remarkably, electron, polariton, and hole systems have 
identical behavior for spatially uniform spin and charge 
densities. Their spin relaxation rates and precession are 
exactly the same. Moreover, there is neither coupling 
between linear and circular (pseudo)spin, nor between 
number/charge density and spin. 

The behavior of spatially inhomogeneous systems is 
more intriguing. An initially unpolarized and non- 
uniform pump will spontaneously produce both circular 
and linear polarization. The resulting circular polariza¬ 
tion has the same magnitude for all three types of par¬ 
ticles, but for holes and polaritons in the EY regime it 
can be zeroed by tuning the magnetic field to a resonance 
condition where the magnetic field matches the spin-orbit 
coupling. The resonance condition differs for holes and 
polaritons. In addition to the circular polarization, the 
linear polarization also can be zeroed in the EY regime 
by tuning the magnetic field to a resonance with the spin- 
orbit coupling, but only in holes. 

We have also seen that an initially polarized pump will 
undergo out-of-plane precession, converting linear polar¬ 
ization to circular polarization and vice versa. Both this 
process and the production of linear polarization from an 
unpolarized pump are strongest for electrons and weak¬ 
est for holes, and produce distinctive dipole, quadrupole, 
and sextupole patterns for the three types of particles. In 
the case of out-of-plane precession, this pattern’s magni¬ 
tude and angular orientation are both sensitive to the 
dimensionless energy splitting C. 

The analytical results presented here were obtained 
with a spatial Fourier transform, i.e. in momentum 
space, and they determine the matrix D~^ which con¬ 
trols time evolution. We now turn to numerical calcu¬ 
lation of the response function $ in real space. These 
numerical results will give us a clear picture of both the 
spatial distributions that are produced from a spatially 
localized pump, and of the spin polarization signals that 
occur before and after the spin relaxation time. 


V. NUMERICAL RESULTS 


In the previous analytical section we calculated the 
diffusion operator D~^{q) in momentum space, which we 
presented in equations [25pT The response function 
in momentum space can not be obtained without numer¬ 
ical calculation, either inverting zw + to obtain 

the frequency response $(a;, or exponentiating 
to obtain the time response Calculation of the 

real-space response functions ^{uj,x) requires a 

further Fourier transform; we discretize on a spatial lat¬ 
tice and use fast Fourier transforms. The lattice spacing 
a adds an additional smearing of the response function 
$ over an area a^, in addition to our diffusive approxi¬ 
mation which removes ballistic physics at scales smaller 
than the scattering length 1. If a < I this lattice-induced 




FIG. 2: (Color online.) Polarization Precession. The black 
lines show the polariton density, and the orange, red, and 
blue lines show linear Sx, Sy and circular Sz polarization. 
Polarization density S is shown in panel a, and polarization 
degree s is shown in panel b. The initial polaritons are 100% 
linearly x polarized. In panel a the initial polariton density is 
spatially uniform producing strictly in-plane precession, while 
in panel b it has wave-vector ql = l/vTO which allows out-of- 
plane precession. Panel a shows in-plane spin precession - out- 
of-phase oscillations of the Sx, Sy polarizations. Solid lines 
show precession with both Zeeman and spin-orbit terms {9b = 
Stt/S), and dashed lines show precession with only a Zeeman 
term {9b = 0.) Panel b shows out-of-plane precession, which 
produces circular Sz polarization. = 0.1, x = h/EpT = 
l/dzr, Tdecay = 50t, and r = 4ps. In panel h 9 b = Stt/S. 


smearing is insignificant. We employ periodic bound¬ 
ary conditions for numerical convenience, combined with 
large sample sizes. In experimental realizations if po¬ 
laritons reach the device edges before they decay, then 
they escape from the edges of the device. The behav¬ 
ior of electrons and holes at device edges is sensitive to 
attached leads. In the diffusive regime these effects are 
not significant for the large sample sizes which we employ 
here. 


A. Basic Spin Precession and Decay 

We begin our numerical results with two graphs that 
illustrate quantitatively basic features of the diffusion 
equations in momentum space. First, Figure stud¬ 
ies the spin precession and decay of systems with fixed 
wave-vector. Panel[^ illustrates precession of an initially 
x-polarized density which is spatially uniform (zero mo¬ 
mentum g = 0) and therefore shows only in-plane and 
no out-of-plane precession. The in-plane precession is 
caused by the Zeeman field oriented perpendicularly to 
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FIG. 3: (Color online.) Polarization (Spin) Decay Times 
Tp as a function of the dimensionless spin splitting ( - both 
for circular polarization (unbroken lines) and for linear 
Sx polarization (dashed lines.) We plot Tp at four values of 
9b ~ 0,7r/6,7r/3,7r/2. In the EY ^ 1/2 regime the decay 
time is constant, and in the DP 1/2 regime it scales with 
C“^. X = h/EpT = I/Itt, Tdecay = 200ps = 50r, and r = Aps. 


the sample, and is manifested as out of phase oscillations 
of the Sx and Sy linear polarizations (orange and red 
lines). Dashed lines show precession with only a Zeeman 
splitting, and solid lines show precession with both a Zee- 
man splitting and a spin-orbit term. The black solid line 
shows the overall polariton density, which decays expo¬ 
nentially due to the effects of the finite polariton lifetime. 
Identical results are obtained for electrons and for holes, 
which at uniform density have identical spin relaxation 
dynamics. Of course, because the lifetime of electron and 
hole charge is infinite, their total density is conserved. 

Panel b of Figure shows out-of-plane spin preces¬ 
sion, where circular polarization is generated from 
the initial x polarized pump. Unlike panel a where we 
plotted the polarization densities Sx, Sy, here we plot 
the polarization degrees Sx,Sy,Sz, which are the polar¬ 
ization densities divided by the polariton density I. We 
will plot this observable s in all subsequent graphs be¬ 
cause it shows out-of-plane precession more clearly than 
the polarization densities S and moreover can be mea¬ 
sured easily in experiments. In order to obtain out-of- 
plane precession, the polariton distribution must have 
non-zero wave-vector, which we have set io q = x/ZvTO. 
The largest polarization degrees are found at large times 
t T. At time t = 50r the polariton density is about 
3.0% of the initial value, the initial 100% x polarization 
degree has reversed sign and is 41%, and a 13% circular 
polarization degree has been generated. 

In Figure]^ we show the polarization decay times of a 
spatially uniform distribution as a function of (/. Identical 
results are obtained for electrons, polaritons, and holes, 
which in the case of spatially uniform systems have iden¬ 
tical spin relaxation dynamics. We show Tp at four values 
oi 9 b- The D’yakonov-Perel’ ^ 1/2 regime is shown 

on the left side, and the Elliott-Yafet regime C 3> 1/2 
regime is shown on the right side. In the EY regime Tp 
is locked to the scattering time r, and in the DP regime 



r / '/Dt r / y/DT 


FIG. 4: (Golor online.) Radial Distribution of the Polari¬ 
ton Density and Polarization Degree. The left panels show 
the time response pit, x) at time t to a transient pulse, while 
the right panels show the frequency response to an oscillat¬ 
ing signal with period T. All panels show the signal’s radial 
dependence, after performing an rms average over the po¬ 
lar angle />. The radial coordinate is rescaled by the diffusion 
length ro = y/Dt. The upper panels plot the polariton density 
(in arbitrary units) produced by an unpolarized pump, 
the middle panels plot the circular polarization degree (s^),^ 
produced by out-of-plane precession from a linearly polarized 
pump, and the bottom panels plot the linear polarization de¬ 
gree {sx) 4 , produced by an unpolarized pump. For polaritons 
Tdecay = 50r. ( = 0.1, 6b = Stt/S, X = ft/EpT = l/dTT, T = 
Aps. 

it scales with At very small C the polariton decay 
time Tdecay = 50t caps Tp. In electron and hole systems 
this cap is removed and Tp continues to scale with 
even at very small C- 

The straight solid blue line in Eigurej^shows that when 
9b = 0, i.e. when the spin-orbit coupling is zero, the 
circular Sz polarization lifetime is equal to the polariton 
lifetime Tdecay ■ In this case Uz commutes with the Hamil¬ 
tonian and circular polarization is conserved. Therefore 
we find large accumulations of circular Sz polarization 
when 9b is small, i.e. when the spin-orbit coupling is 
small. 


B. The Polarization’s Dependence on Radius 

Next we begin our study of the response to a spatially 
localized pump, which requires numerical evaluation of 
equations and In order to learn about typical 
length scales. Figure ^ focuses on the radial dependence 
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of the polarization. We present radial data for polaritons 
only because the electron and hole results are quantita¬ 
tively similar and qualitatively identical. The similarity 
of the radial dependence between the three particles, and 
indeed all of the qualitative and quantitative results seen 
in Figure]^ and discussed here, was inaccessible without 
numerical calculations. 

The left panels of Figure]^ show the temporal response 
at times t = 5r, 20t, 80t after the pulse, and the right 
panels show the frequency response to an oscillating sig¬ 
nal with periods T = 5t, 20r, 80r. The blue lines in 
the right panels show the response to a static pump. 
In the t,T = 5t, 20t calculations we used a system of 
size L = W = 51.21 and a lattice spacing a = 0.21. 
In the t,T = 80r calculations L = W = 102.4 Z and 
a = 0.4:1, and in the static calculation L = W = 204.8 Z 
and a = 0.41. Since a I, lattice-spacing effects are 
insignificant. We rms-averaged the intensity I and polar¬ 
izations Sx,Sy,Sz over the polar angle (j). The averaged 
polarization degrees {Sx, Sy, 8^)4, are calculated by divid¬ 
ing the rms-average of Sx,Sy,Sz by the rms-average of 
/. 


The upper panels of Figure]^ plot the polariton den¬ 
sity I, which is independent of and also independent 
oi 9 b and of the starting polarization, except for small 
variations when 9b = 0. In panel a the time-response 
curves are simple quadratics, and they match perfectly 
after rescaling r by y/Dt, which shows that the polariton 
density is simply an expanding Gaussian. Similarly in 
panel b the frequency response curves coincide well after 
rescaling by rg = V DT. The frequency response decays 
exponentially with r for all r/rg > 0.3. (We obtain a 
good match for the static data by rescaling with the ar¬ 
tificial value T = 1280ps. This value of T may be tied 
to the decay time Tdecay = 200ps.) The density’s almost 
perfect insensitivity to (^,9b, and time t, T, after rescal¬ 
ing the radius r, supports our use of intensity variables 
Sx, Sy, Sz in the remaining graphs. Because / has a sim¬ 
ple profile, dividing by I does not obscure or complicate 
the physics revealed in the following graphs. 

The density’s form is very similar to that of a sim¬ 
ple particle density undergoing diffusion without spin- 
orbit coupling or spin precession. In this case the dif¬ 
fusion operator is just Drq^, where D is the diffusion 
constant, and the Fourier transform to position space 
can be performed analytically. For the time response 
4>{t,x) (Equation 21) one finds an expanding Gaussian 
{2TTDt)~^ exp{—r'^/2Dt). For the frequency response 
$(a;,a:) (Equation 22) one finds a modified Bessel func¬ 
tion 2TT\Ko{ry/uJ/D)\, which decays exponentially with 
r if 0.2 < r^JujjD. Our numerical results are very similar 
to these analytical forms. 

The middle panels c and d show out-of-plane preces¬ 
sion from a linearly x polarized pump to circular Sz po¬ 
larization degree. The static response (blue line) shown 
in panel d displays spatial oscillations which are caused 
by the regularization that we imposed on the diffusion 
operator at large wave-vectors ql > 1/^[2. However we 


have checked that the overall shape and magnitude of 
the static response is insensitive to the cutoff. In addi¬ 
tion, our calculations using other cutoffs show that the 
frequency response at T = 5t, 20r, 80t is sensitive to the 
cutoff at small radii r/V DT < 2, and other cutoffs not 
shown here can show substantially larger frequency re¬ 
sponses at small r. Outside of this window the frequency 
response is insensitive to the cutoff. The temporal re¬ 
sponse is insensitive to the cutoff everywhere. 

Panel c shows that at t = 5r, 20t conversion rates 
can reach and exceed 10“^. Panel d shows that similar 
numbers are obtained from the response to oscillating 
and static pumps. These numbers decrease substantially 
at larger times or periods t,T = 80t, and can be tuned 
by varying C and 9b. An increased polarization degree 
can also be obtained by moving to larger radii r, at the 
cost of reducing the polariton density. 

Panels e and f show the linear Sx polarization degree 
produced in response to an unpolarized pump. These 
curves are very similar to the data in panels c and d, 
but numerically they are about an order of magnitude 
smaller. This difference in magnitude is caused by the 
factor of X, the ratio of the energy splitting to the pump 
energy Ep, which controls angular variations in the mo¬ 
mentum kp. X is l/47r « 0.08 in our calculations, x 
occurs once in the coupling kqi of polariton density to 
circular polarization, signalling that this process requires 
appreciable angular variations oi kp. x does not occur at 
all in the coupling of linear to circular polarization kic] 
therefore panels e and f are smaller by a factor of x- 

Our time response and frequency response results show 
that the polarization degree - the ratio of the polariza¬ 
tion signal to the density - generally increases with r/rg. 
This is a universal trend visible in all of our data, for 
electrons, polaritons, and holes, and for all values of the 
model parameters, and is seen in panes c,d,e, and f. To 
be clear, the polarization signal itself decreases with r. 
However its rate of decrease is smaller than that of the 
number/charge density N, which results in the observed 
increase in r seen throughout Figure At large enough 
r the polarization degree invariably exceeds 1, as seen for 
example in pane d. However this unphysical result occurs 
only when the signal intensity is smaller than the peak 
value by a factor of 10““^ or less; it is not physically mea¬ 
surable. This unphysical result is a mathematical artifact 
of our diffusive approximation, which does not correctly 
capture the tails of probability distributions. In the re¬ 
maining graphs we restrict r to values where the diffusive 
approximation remains valid. 


C. Angular structure in position space. 

We turn to the polarization’s angular distribution in 
real space, which naturally complements the previous 
section’s analysis of the radial dependence. Our numer¬ 
ical results will confirm the analytical results about the 
multipole patterns of electrons, polaritons, and holes, and 
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FIG. 5: (Color online.) Sz circular polarization degree in real 
space, at two values of (, and two values of time t. The left 
panels and right panels are t = 5 t and t = 20t respectively, 
and the upper and lower panels are at = 0.25 and ^ = 
0.5 respectively. The pump is linearly x polarized. As ^ is 
increased the polarization’s quadrupole distribution rotates. 
As time increases the quadrupole pattern rotates and deforms 
into a spiral. 9b = 57r/8, x = h/EpX = l/lrr, Tdecay = 50r, 
and T = 4ps. 


FIG. 6: (Color online.) Sz spin polarization degree in real 
space, for electrons (upper panels) and holes (lower panels). 
The left panels and right panels are t = 5t and t = 20t 
respectively. The pump is linearly x polarized. Electrons 
show a dipole pattern, and holes show a sextupole pattern. 
As time increases the pattern rotates and deforms into a spiral 
in response to the magnetic field. ( = 0.25,03 = 57r/8, x = 
h/EpT = l/47r, Tdecay = oo, and T = 4ps. 


their rotation with Furthermore, we will show that the 
multipoles twist into spirals as time progresses, and we 
will give quantitative results about the signal strengths. 

Figure]^ shows the pattern of circular polarization de¬ 
gree Sz = Sz/I which out-of-plane precession produces 
starting from a linearly polarized pump. Our analytical 
results indicate that direct conversion of linear to circular 
polarization in the EY regime requires both a spin-orbit 
term and a Zeeman splitting, so we make both terms 
strong with 9b = Stt/S. Here and in the other radial 
plots (Figures 1^ and 10) the lattice spacing is 0.2/ and 
the system size is 51.2/. Previous theoretical and exper¬ 
imental articles in the ballistic regime have shown that 
out-of-plane precession produces a quadrupole pattern 
with two opposite lobes of positive circular polarization, 
two negative lobes, and the Langbein cross of zero polar¬ 
ization dividing these lobes.^^HH Further work, always in 
the ballistic regime, showed that when a magnetic field 
is introduced the quadrupole pattern twists into a spi¬ 
ral pattern^ with its orientation becoming a function of 
radius This rotation occurs because in the ballistic 
regime the radial position r is proportional to time t, and 
the external Zeeman field causes in-plane precession at a 
fixed frequency, which in ballistic fight maps to a preces¬ 
sion length. The combination of the quadrupole pattern 
with spin in-plane precession causes the spiral. 

As we discussed earlier, our analytical results for the 
linear-circular coupling k/c do predict a quadrupole pat¬ 
tern, and in addition indicate that the quadrupole pat¬ 
tern rotates by 7r/4 as the ( parameter is shifted from the 
DP regime to the EY regime. The left panels of Eigurej^ 


show the circular polarization degree at time t = 5t, at 
C = 0.25 (upper left) and at ^ = 0.5 (lower left.) These 
panels show a clear quadrupole pattern, and confirm that 
the pattern rotates as C is changed. Similarly to the bal¬ 
listic regime, the rotation is caused by in-plane precession 
superimposed on the multilobe pattern. However - this 
is a hallmark which can be used to distinguish the dif¬ 
fusive regime from the ballistic regime - no concentric 
circles are observed in the diffusive regime. Surprisingly, 
even though there is a Zeeman term the polarization pat¬ 
tern is not twisted into a spiral, unlike previous results 
in the ballistic regime at t < r. It seems that the spi¬ 
ral effect does not occur at this time scale. However the 
right panels of Figure at t = 20t do show clear spiral 
patterns at both ( = 0.25 (upper right) and at C = 0-5 
(lower right.) Appendixshows that these patterns are 
mildly sensitive to the cutoff but the main features are 
cutoff-insensitive. Panel d with its very small polariza¬ 
tion degree is the exception: its pattern does depend on 
the cutoff. In summary, the circular polarization pattern 
deforms from a simple quadrupole pattern to a spiral as 
time progresses, and its angular orientation is sensitive 
to both C and time t. The most surprising aspect of these 
results is that no spiral pattern is seen at time t = 5r, 
unlike ballistic results. Experimental observation of the 
spiral pattern’s temporal evolution from a simple multi¬ 
pole to a spiral may be a useful means of measuring the 
scattering time in individual devices. 

Figure shows the angular patterns produced by elec¬ 
trons (top panels) and holes (lower panels), which have 
been little discussed in the previous literature. As ex- 
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pected, electrons display a dipole pattern in real space, 
and holes show a sextupole pattern. Comparison of the 
left t = 5t and right t = 20r panels shows again a time 
dependence - the initial patterns twist into spirals in re¬ 
sponse to the magnetic field. In results not shown here, 
we also have confirmed that the angular orientation of 
the electron dipole pattern and the hole sextupole pat¬ 
tern rotate in response to variations in C, similarly to the 
polariton’s C dependence. Probably the most interest¬ 
ing aspect of these results is that the observed real-space 
spin pattern (dipole, quadrupole, sextupole) is a simple 
signature that can be used to determine the dominant 
type of spin-orbit interaction. Moreover, it is likely that 
in materials with both linear and cubic spin-orbit terms, 
a numerical analysis of experimental spatial spin distri¬ 
butions could distinguish both terms in the spin-orbit 
interaction. 


D. The Polarization’s Temporal Dependence 

In the previous sections we analyzed the spin’s spatial 
profile; both its radial and angular dependence. Here we 
turn to the time dependence of the spin polarization de¬ 
gree, which is determined by the response function <I‘(t), 
the exponential of the diffusion operator. In Figure 
we examine the time dependence of the spin polarization 
degree of electrons (top panels), polaritons (middle), and 
holes (lower panels), both in the DP regime (left panels) 
and in the EY regime (right panels.) In each panel the 
dashed and solid lines show respectively the linear Sx po¬ 
larization degree and the circular Sz polarization degree. 
Blue lines show the Sx , Sz polarization degrees produced 
by an unpolarized pump, and orange lines show the po¬ 
larization degrees caused by out-of-plane precession from 
a polarized Dump. The lattice spacing a, both here and 
in Figures 1^ and is a = 0.2 ro = 0.2-i/Dt. Since the 
lattice spacing is effectively the size of the spatially lo¬ 
calized source pump, this means that we are varying the 
pump size, with maximum width at t = 125 r of about 
a « 1.6 1. However the effect of this smearing remains 
small in our plots, which show the spin polarization de¬ 
gree at much large distances r = Stq. The sample width 
and length are 12.8ro, large enough to make finite size 
effects quite small. We set 0 b = Stt/S, i.e. both the Zee- 
man splitting and the spin-orbit splitting are significant, 
but not precisely in balance. This allows both in-plane 
and out-of-plane precession to occur, and therefore gives 
larger signal strengths than those seen at, for instance, 
9b = 0 or 9 b = '!^I2.. At times up t = 4 t the data may 
be influenced by ballistic effects which we omitted in our 
diffusive approximation, and in some cases is sensitive to 
the cutoff. In particular, with other cutoffs some spin sig¬ 
nals are increased. These ballistic/cutoff effects at small 
times have little effect on the data after t = At. 

The spin signals in Figurej^generally peak at t = 4—8r 
in the DP regime and at t = 1.75 — 6r in the EY regime. 
Obviously this time scale is determined by the scatter- 



FIG. 7: (Color online.) Temporal dependence in electrons 
(top), polaritons (middle), and holes (bottom). ( = 0.1 in 
the left panels and ^ = 4 in the right panels. Blue lines show 
the polarization degree produced by an unpolarized pump, 
and orange lines show the polarization degree caused by out- 
of-plane precession from polarized pump. Dashed lines show 
linear Sx polarization degree, and solid lines show circular Sz 
polarization degree. The data is rms-averaged over the polar 
angle cj>. For polaritons Tdecay = 50r and for electrons and 
holes Tdecay = oo. Tj\fDt is kept fixed at 3, 9b = Stt/S, x = 
IjEpT = l/47r, r = Aps. 


ing time t. At this time scale, in the DP regime, the 
out-of-plane precession from an initially polarized pump 
(orange lines) is dominant, an order of magnitude larger 
than polarization generated from an unpolarized pump 
(blue lines). In the EY regime the two processes have 
similar magnitudes, because any initial polarization de¬ 
cays extremely quickly. The peak polarization degree is 
greatest for electrons: 61% in the DP regime, and 32% 
in the EY regime. Eor polaritons it is 43% and 8% re¬ 
spectively for the DP and EY regimes, and for holes it is 
23% and 11%. These peak values are independent of the 
cutoff in the DP regime. They are more sensitive in the 
EY regime where the peak occurs earlier. As a general 
trend, the polarization signals are invariably larger in the 
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DP regime than in the EY regime. 

Starting at t = 4 — 8r in the DP regime and at 
t = 1.75 — 6r in the EY regime, the polarization degree 
decreases uniformly. At the same time the dashed lines 
merge with each other, and solid lines also merge with 
each other. This means that the system loses its mem¬ 
ory of the initial polarization, and is sensitive only to the 
original polariton intensity. The memory loss implies a 
reversal in the dominant mechanism of spin production: 
before the memory loss out-of-plane precession is domi¬ 
nant, while after the memory loss it is replaced by spin 
generation from an unpolarized source. 

In the EY regime (right panels) the memory loss com¬ 
pletes before t = lOr, because the spin is tightly locked to 
the momentum and is randomized at every scattering. In 
the DP regime (left panels) the memory loss lasts much 
longer - lOOr or longer - because the spin is little affected 
by scattering and instead is lost by spin-orbit precession 
around a randomized precession axis. This is consistent 
with the analytical expectation that in the DP regime 
the spin relaxation time scales with t/^^. In our case 
with ^ = 0.1 the memory loss occurs a factor of 10 later 
than the time of the peak in the polarization degree. 

Surprisingly, after memory of the initial polarization 
is completely lost all polarization degree signals do not 
decay exponentially, but instead as powers of the inverse 
time (straight lines on our log-log plots). This power law 
decay allows the polarization degree to persist to quite 
long times: at t = lOOr the Sz polarization degree for 
polaritons is 0.5% in the DP regime, and for electrons it 
is 7%. In the EY regime the t — lOOr polarization 
degree is smaller, but is still 2% in the case of electrons. 

It is also interesting that after the memory loss the 
polarization degree is largest for electrons and smallest 
for holes, because kqi’s out-of-plane precession between 
Sx, Sy and Sz is linear for electrons, quadratic for po¬ 
laritons, and cubic for holes. In particular, in polaritons 
the Sx and Sz signals are almost the same because in 
this case all couplings are quadratic in the momentum. 
In contrast, the perpendicular Sz polarization degree is 
roughly the same for all particles, in keeping with the doz 
coupling which is the same up to a numerical constant. 


E. Dependence on the Dimensionless Energy 
Splitting 

In Eigure|^we examine the polarization degree’s de¬ 
pendence on the energy splitting which is the ratio 
of the scattering length I to the precession length Iprcsn- 
We will confirm the analytical results’ linear dependence 
on C both in the DP and EY regimes. Going beyond 
the analytical results, we will reveal the behavior at the 
transition between DP and EY regimes, give quantita¬ 
tive magnitudes of the signals, and find an interesting 
feature of polaritons in the EY regime. Again we have 
set 9b = 57r/8 to allow both in-plane and out-of-plane 
precession to occur, and to give larger signal strengths. 



FIG. 8: (Color online.) Dependence on the dimensionless 
energy splitting ( = Ip/lprcsn. in electrons (top), polaritons 
(middle), and holes (bottom). Similarly to Figure blue 
lines show the polarization degree produced by an unpolar¬ 
ized pump, and orange lines show the polarization degree 
caused by out-of-plane precession from a polarized pump. 
Dashed lines show linear Sx polarization degree, and solid 
lines show circular Sz polarization degree. The data is rms- 
averaged over the polar angle <j)- Here r/'/Dt is kept fixed at 
2>,t = 20r, 6b = 57r/8, y = 'VjEpT — l/dvr, r = Aps. 


The time is fixed at t = 20r, long after the first scatter¬ 
ing. We plot the same observables seen in Figure with 
the same parameter values. The data between ^ = 0.25 
and C = 1 is mildly dependent on the cutoff, but the qual¬ 
itative trends are insensitive to the cutoff and to ballistic 
physics. 

The polarization degree signals show a transition which 
begins near C = 0.1 and ends near C = 1/2, i.e. at the 
transition from the DP regime to the EY regime. This 
transition is controlled by the onset of spin relaxation, 
which occurs at t oc t/(/^ in the DP regime, suggest- 
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ing that the transition in C scales with . At small 
C < 0.1, i.e. deep in the D’yakonov-Perel’ regime, the 
out-of-plane precession signals (orange lines) are largest 
and tend to have very similar values. These signals rise 
steadily with C, in agreement with the linear-circular cou¬ 
pling Kjc which is linear in C in the DP regime. Produc¬ 
tion of linear polarization and circular polarization 
from an unpolarized pump (the blue lines) is smaller. In 
polaritons and holes this process is linear in C, again in 
agreement with the kqi and doz couplings which are linear 
in C. 

Near ( = 0.1 the out-of-plane precession signals (or¬ 
ange lines) reach peak values of 43% in electrons, 24% in 
polaritons, and 8% in holes. Between C = 0.1 and C = 1/2 
the out-of-plane (orange) signals drop precipitously, and 
for polaritons and holes the signals from an unpolarized 
pump (blue lines) also decrease. Near the transition to 
the EY regime at ^ = 0.5 the solid lines merge with each 
other, as do also the dashed lines. The mergers indicate 
that the system loses its memory of the pump polariza¬ 
tion; it is sensitive only to the initial pump intensity, not 
to the initial polarization. This is because at < = 20t a 
system in the EY regime is insensitive to its starting po¬ 
larization. In the EY regime all polarization signals are 
linear in C, as expected from our couplings kqp doz, 
which are all linear in in the EY regime. We again see 
that the Sx polarization degree is largest for electrons and 
smallest for holes, because koz’s out-of-plane precession 
between Sx, Sy and Sz is linear for electrons, quadratic 
for polaritons, and cubic for holes. In contrast, the per¬ 
pendicular Sz polarization degree is roughly the same for 
electrons, polaritons, and holes. Interestingly, polaritons 
in the EY regime give almost exactly equal magnitudes 
for all four of the signals plotted here. 


F. Dependence on 9 b , the balance between 
spin-orbit and Zeeman terms. 

Figure examines the polarization degree’s depen¬ 
dence on the angle 0 b, which describes the balance be¬ 
tween the spin-orbit interaction and the Zeeman split¬ 
ting. When 9b =0 there is only a Zeeman splitting and 
no spin-orbit interaction, and when 9b = 7r/2 this is re¬ 
versed. This section confirms the 9b dependence which 
we discussed in the analytical section; its main value is 
to illustrate the analytical results and give quantitative 
magnitudes. The time is fixed at t = 20t, long after the 
first scattering, and the data is almost completely insen¬ 
sitive to cutoff effects. We plot the same observables seen 
in Figures and with the same parameter values. 

If one flips the sign of both the spin-orbit term and 
the Zeeman splitting, or equivalently adds tt to 9b, 
then the polarization signal may change its overall sign 
but its magnitude must be unchanged. Therefore our 
plots of the rms-averaged signal are symmetric under 
9b —i' 9b + TT, and we plot only the interval 9b = [0,7r]. 
In addition we find that three of the four signals in Fig- 



FIG. 9: (Color online.) Dependence on 9b in electrons (top), 
polaritons (middle), and holes (bottom). = 0.1 in the left 
panels and (^ = 4 in the right panels. Blue lines show the 
polarization degree produced by an unpolarized pump, and 
orange lines show the polarization degree caused by out-of¬ 
plane precession from a polarized pump. Dashed lines show 
linear Sx polarization degree, and solid lines show circular 
polarization degree. The data is rms-averaged over the polar 
angle 4>- rjy/Dt is kept fixed at 3, t = 80ps = 20r, x = 
hjEpT = l/47r, T = 4ps. 


urej^are symmetric under flips around 9b = tt 12. These 
signals are independent of the sign of the Zeeman term, 
which distinguishes between clockwise and counterclock¬ 
wise precession in the device plane. The one exception 
occurs in the EY regime, where the linearly polarized Sx 
response to a circularly polarized pump (dashed or¬ 
ange lines) is not symmetric under flips about tt/2 , i.e. it 
is sensitive to the Zeeman term’s sign, as a result of the 
pump’s orientation parallel to the z axis. This sensitivity 
occurs only if the system retains memory of the original 
polarization, i.e. either at short times, or even at long 
times when 9b is close to zero as seen here. However, 
even this signal is symmetric under changes at the same 
time of both the sign of the pump’s Sz polarization and 
also sign of the Zeeman term. 
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A first examination of Figure shows that the polar¬ 
ization signals are strongly sensitive to 9 b- In many of 
our other graphs we have fixed 9b = Stt/S in order to 
avoid the deep minima that occur at 9b = 0,7r/2, and 
elsewhere. We also see again that in the EY regime (right 
panels) the solid lines coincide, and the dashed lines also 
coincide, except near 9b = 0 ,7^ where is conserved. As 
noted before, this means that the system has forgotten 
its initial state. 

As discussed earlier, a special feature of electrons is 
that in the DP regime the conversion from unpolarized 
pumps to linear polarization scales with a factor of 

smaller than the conversion to circular polarization, 
and also a factor of smaller than the same signals in 
polaritons and holes. In Figure we see that the unpo¬ 
larized to linear signal (the dashed blue line) is smaller, 
but not by the factor of = 10“^ that might be ex¬ 
pected. The reason is that here we are seeing the result 
of a two step process: (a) conversion of to circular polar¬ 
ization, and (b) precession from circular to linear polar¬ 
ization. This also explains why the unpolarized to linear 
signal is zero at 9b = 7r/2, i.e. when there is no magnetic 
field. 

Broadly speaking, the details of Figure confirm fea¬ 
tures of the 9b dependence which we found in our anal¬ 
ysis of the couplings In particular: 

• No linearly polarized signal is produced when there 
is no spin-orbit coupling, i.e. at 9b = 0. This is 
universal for electrons, polaritons, and holes, and 
occurs because a Zeeman field oriented along the 
z axis is unable to create a linear polarization, be- 
cuase it is unable break the rotational symmetry in 
the xy plane. 

• At = 0 the absence of a spin-orbit term ensures 
that the lifetime of circular Sz polarization is infi¬ 
nite. Therefore at 0 b = 0 we find large circular Sz 
polarization degrees. The Sz magnitude at 0 b = 0 
is close to xC = C/47’‘ « 0.08 x Q, and is insensitive 
to C and to any initial linear Sx polarization. It is 
also the same for electrons, polaritons, and holes. 

• In the DP regime (left panels), regardless of the 
value of 0B, a linearly x polarized pulse always gen¬ 
erates Sz polarization. When there is no spin-orbit 
coupling, i.e. 0 b = 0, the Sz polarization is gen¬ 
erated by conversion from the pump’s initial num¬ 
ber/charge density, and its magnitude is of order 
XC- When there is no magnetic field, i.e. 0 b = 7r/2, 
the spin-orbit coupling causes out-of-plane preces¬ 
sion and produces the Sz polarization. In this case 
its magnitude reaches 48% in electrons, 27% in po¬ 
laritons, and 9% in holes. 

The most interesting aspects of Figure are the differ¬ 
ences between electrons, polaritons, and holes: 

• In the EY regime (right panels) polaritons exhibit 
no polarization signal at all unless there is a Zee- 
man term, so at 0 b = 7^/2 all signals are zero. This 


contrasts with electrons and holes in the EY regime, 
where a Zeeman term is not necessary for produc¬ 
ing an Sx signal, and is required only for producing 
Sz- This means that for electrons and holes the 
dashed lines are nonzero at 0 b = 7r/2. 

• In both the EY and the DP regime an unpolar¬ 
ized pump does not produce a linearly polarized Sx 
signal unless a magnetic field is introduced, so at 
9b = 7r/2 the solid blue lines are zero. 

• Setting 0B = 7r/4, 37r/4 creates a resonance condi¬ 
tion where the Zeeman and spin-orbit terms have 
precisely the same magnitude, i.e. in polaritons 
\bz\ = |Afc|.|. For polaritons in the EY regime this 
resonance condition zeroes the generation of circu¬ 
lar Sz polarization f (graphed as solid blue and or¬ 
ange lines), causing sharp dips at 9b = 7r/4, 37r/4. 
Holes in the EY regime also display a resonance 
effect, with a shifted resonance which occurs at 
\bz\ = v^lAfc^j. In addition, holes manifest a sim¬ 
ilar resonance condition in the generation of Sx lin¬ 
ear polarization from an unpolarized pump, shown 
as dashed blue and orange lines. In contrast elec¬ 
trons do not display any resonance at any value of 
bz- The resonances shown here may be useful for 
measurements of the spin-orbit strength in polari- 
ton and hole systems in the EY regime. 

G. Summary of the Numerics 

The main value of our numerical results, as opposed 
to our analytical formulae, is that they offer a clear pic¬ 
ture of the spatial spin patterns in these systems, and 
also of their time evolution. Our spatial data on the 
number/charge density showed that the frequency re¬ 
sponse always follows a simple decaying exponential in 
real space, while the temporal response follows a Gaus¬ 
sian. Moreover, our data showed that changing the time 
parameter, or indeed any other parameter, has no influ¬ 
ence on the number/charge density other than rescaling 
its exponential or Gaussian profile. 

The spin polarization shows much richer behav¬ 
ior. Electrons, polaritons, and holes produce dipole, 
quadrupole, and sextupole patterns respectively, which 
provides a very simple and strong diagnostic tool for vi¬ 
sually determining the dominant spin-orbit term. Our 
numerical results confirm that these patterns are clearly 
visible in the diffusive regime. In real materials the spin- 
orbit interaction may involve more than one term, for 
instance both linear and cubic terms, and in this case 
numerical analysis of experimental data will reveal the 
relative strength of the spin-orbit terms. Our numeri¬ 
cal results also showed that in the presence of a mag¬ 
netic field, the spin distribution evolves with time from 
a simple multipole into a spiral. This too will be useful 
for measuring the scattering time in individual devices. 
Moreover, we saw that in the diffusive regime there is 
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no pattern of concentric circles around the pump, which 
allows an unambiguous determination of whether the sys¬ 
tem is in the diffusive or ballistic regime. 

The numerical data also gave us rich information about 
the time dependence, and told us how large the polariza¬ 
tion degree can become. We can break the time evo¬ 
lution into two phases: before and after the time when 
the system loses memory of its initial polarization. This 
memory loss occurs at t cx t in the DP regime, and 
before t = lOr in the EY regime. Before memory loss 
occurs, the dominant source of spin dynamics is out-of- 
plane precession which converts Sx,Sy spin to spin, 
and vice versa. With an optimal value of the dimension¬ 
less spin splitting ^ = 0.1, this process yields maximum 
polarization degrees of 23 — 61% in the the DP regime 
at around t = 4 — 8 r. In the EY regime at ^ = 4 the 
maximum polarization degrees are between 8 and 32%, 
occur at t = 1.75 — 6r, and grow linearly with These 
maximum values are largest for electrons. After the max¬ 
imum, spin relaxation causes a steady decrease in the 
polarization degree. 

As the system loses its memory of the initial polar¬ 
ization, the dominant mechanism producing spin polar¬ 
ization changes from spin precession to direct conversion 
from the number/charge density. In this phase the spin 
polarization degree decreases as a power of the time, and 
in particular the linear Sx polarization degree is largest 
for electrons and weakest for holes. The circular polar¬ 
ization degree is about the same for all three particles. 
Even at t = IOOt the polarization degree for electrons 
can remain as high as 7% in the DP regime, or 2% in the 
EY regime. 

Lastly, our numerical results about the C dependence 
gave the transition between DP and EY regimes (not 
shown in the analytical results), and gave an optimal 
value of C- Our 9b results largely illustrated and con¬ 
firmed the analytical results on this parameter. 


VI. DISCUSSION AND CONCLUSIONS 

After proper rescaling of time and distance, in the 
diffusive regime, spatially homogeneous distributions of 
electrons, polaritons, and holes act identically. They 
all have the same spin relaxation times and the same 
spin precession rates, and in all cases their number 
charge/density shows the same spatial profile spreading 
away from the pump. This is perhaps disappointing, 
since one might wish to see at this level some clear dif¬ 
ference between the bosonic polariton and the fermionic 
electron and hole, similarly to the difference between 
Bose-Einstein condensation and fermionic repulsion that 
is seen at large particle densities. At small particle densi¬ 
ties, in the non-interacting limit, the perturbative treat¬ 
ment utilized here shows no such signal. There is still 
some possibility that a Berry phase or other quantum 
effect could distinguish between the three particles, but 
there is no such effect even in the weak localization signal. 


which is determined by relaxation times that are mathe¬ 
matically identical to the spin relaxation times calculated 
here. Any such quantum effect will be rather subtle. 

However, using a spatially localized pump will imme¬ 
diately reveal very large differences between the three 
particles. The most important difference is their multi¬ 
pole spatial pattern, which has two lobes for electrons, 
four for polaritons, and six for holes. This spatial distri¬ 
bution provides an unambiguous signal to experimental¬ 
ists for determining immediately which spin-orbit term 
is dominant, and also provides the raw data for quanti¬ 
tative analysis of the strengths of competing spin-orbit 
terms. It is a very strong signal: if the pump is initially 
linearly polarized, then in the DP regime a multipole pat¬ 
tern in the circular polarization can be observed, with a 
polarization degree of up to 23 — 61%, and can persist 
for long times. In the EY regime the same process can 
produce between 8 and 32% polarization degree. Pro¬ 
duction of linear polarization from an unpolarized pump 
also produces appreciable multipole patterns. 

The three particles also differ in their response to a 
perpendicular magnetic field. For electrons and holes, in 
the EY regime an unpolarized pump always produces a 
multipole pattern of linear polarization, regardless of the 
magnetic field strength. In contrast for polaritons in the 
EY regime no linear polarization is generated from an un¬ 
polarized pump unless there is a magnetic field. More¬ 
over in the EY regime polaritons and holes manifest a 
resonance, where no circular Sz polarization is generated 
from an unpolarized pump, at a specific field strength 
that matches the spin-orbit strength. The resonance oc¬ 
curs at numerically different field strengths for different 
particles, and varies with spin-orbit strength. This will 
allow experimentalists both to measure the spin-orbit 
strength and to confirm again whether the particle is an 
electron, polariton, or hole. 

The last major difference between the three particles is 
in the magnitude of multipole pattern that is produced. 
In electrons, polaritons, and holes this scales respectively 
with the first, second, and third power of the scattering 
length. In the diffusive regime studied here the scatter¬ 
ing length is small compared to the experimental length 
scale, so the multipole is strongest for electrons and weak¬ 
est for holes. One interesting way of testing this trend 
is to compare, at long times, the linear and circular po¬ 
larization produced from an unpolarized pump. The cir¬ 
cular signal is rotationally symmetric and its magnitude 
is about the same for all three particles. In the case of 
electrons the circular signal should be smaller than the 
linearly polarized signal, in the case of holes this is re¬ 
versed, and in the case of polaritons the two signals have 
very similar magnitudes. 
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Appendix A: The diffuson in the general case, for any value of and for linear, quadratic, and cubic 

spin-orbit terms. 


In this article we report on spin dynamics for the following Hamiltonian: 


H — --h hz<Jz Ak^ 

2m 


^iN4n. 
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2m 


+ bzCTz -f A 


0 {kx - iky) 

[kx + iky)^ 0 


N' 


(Al) 


We consider the cases of linear N = 1, quadratic N = 2, and cubic N = 3 spin-orbit interactions. We assume that the 
impurities conserve the (pseudo)-spin quantum number; if this constraint is relaxed then our results on the matrix 
structure of the diffuson will change substantially. 

For the Hamiltonian H, we obtain the diffuson D~^: 
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The constant contributions to the matrix t jTxx-^'cI tzz-,t jTxy^ which determine the evolution of a spatially uniform 
distribution, are independent of A. The diffusion constant dzz is also independent of A. 

As described earlier, S = sin{9B) = Akp / ^ (Akp)"^ + and C = cos{9b) = bz/(Akp)"^ + describe the 
relative strength of the spin-orbit coupling Akp and the Zeeman term bz- 


Appendix B: Numerical Results on the Cutoff Dependence 

As discussed in the text, some of our numerical results are sensitive to physics at length scales shorter than the 
scattering length I, i.e. ballistic physics. These numerical results are changed by our perturbative expansion of the 
response function in powers of ql. We have tested our results by calculating them in three different ways: 
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FIG. 10: (Color online.) Cutoff dependence of our results. Circular Sz spin polarization degree in real space, for cutoff A 
(upper panels) and cutoff B (lower panels). The left panels and right panels are t = 5 t and t = 20r respectively. The pump is 
linearly x polarized. ( = 0.25, 9b = Stt/S, x = h/EpT = l/drr, and r = 4ps. 


• Cutoff A: At small wave-number 0.5 > we use the formulas presented in Appendix]^ At large wave-number 
(ql)'^ > 2 we set kqi, kic, and Dq^ to zero, and we set A = {ql^/2. At intermediate wave numbers 2 > [ql)"^ > 0.5 
we perform a linear interpolation of Kq/, k/o ^ between these two cases. 


Cutoff B: At small wave-number 0.5 > {qVf' we use the formulas presented in Appendix 


{qly > 2 we use the expressions for the Dyakonov-Perel limit, as listed in equations 2^ 27 28 


i^ At 


At large w ave- num ber 
At 


29 and 31 


intermediate wave numbers 2 > (gZ)^ > 0.5 we perform a linear interpolation between these two cases. 


Cutoff C: We use only the formulas presented in Appendix without any changes, for all values of ql. 


Since Cutoff C leaves the analytical formulas unchanged, it retains the original instability near C, = 1/2. Both Cutoffs 
A and B do not show instability for electrons and polaritons. Holes are stable only with Cutoff A and are unstable 
with both Cutoffs B and C, presumably because our analytical expressions for holes include cubic {ql)^ terms which 
are more unstable than the quadratic terms seen in electrons and polaritons. 

Figure [T0| illustrates the differences between two cutoffs. The upper panels are obtained using Cutoff A, and the 
lower panels with Cutoff B. Very similar patterns and magnitudes are seen, with some mild distortions. We have 
performed similar comparisons, with all three cutoffs, for all of our real-space data. In cases where different cutoffs 
give different results, we have discussed this in the main text. 


Appendix C: Derivation of the kq; and doz Couplings 


Here we report the details of our derivation of the coupling kqi between charge and in-plane Sx, Sy spin, and also in 
the coupling doz between charge and out-of-plane spin. We have performed the calculation systematically in a very 
general mathematica script which allows calculation to quite high orders in the various small parameters. However 
here we do not discuss the script, and instead show in analytical formulas how the derivation develops. We compute 
the doz coupling which is equal to —/o3j and the k^i coupling which is equal —/qi and —/02, where lij is the scattering 


operator. We begin with equation 24 which writes lij as the integral 


dTTZZT 


dkTr{ G^ik — q/2, Ep)aiG^{k + q/2, Ep) aj) 


(Cl) 


Here G^ and G^ are the disorder-averaged single-particle Green’s functions, q is the diffuson wave-vector, and is 
the density of states. The trace is taken over the spin indices of G'^, G^, ai, and aj, which are all 2 x 2 matrices in 
(pseudo)spin space. The advanced and retarded Green’s functions can be resolved into separate contributions from 
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the two spin states s = ±1, which have energies E{s, k): 
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Performing the trace and applying the identity (ab)~^ = (a~^ — b~^)(b — a)~^ quickly obtains 
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(C3) 


In the last line we replaced {Ep — E{s, k — q/2) — ih/2T) ^ by its imaginary part, which is taken to be the Dirac delta 
function nx5{Ep — E{s,k — q/2)). This standard approximation is justified if the scattering time is long compared 
to the kinetic time scale, i.e. oc x 1. Physically it means that we concern ourselves only with processes that 
occur on the elastic scattering circle for polaritons, or on the Fermi surface for electrons and holes. After these steps 
the scattering time r occurs only in the denominator of Y. Powers of dimensionless energy splitting C = EprcsnT/h, 
because they contain r, can come only from the denominator of Y. 

Multiplying {sXi{k — q/2) + s Xi{k + q/2)) by the two terms in Y and making appropriate shifts fc —>■ fc± q/2 
obtains 
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At this point we expand the Xi functions in powers of q. Each power of q is implicitly accompanied by a power of 
X/kp., and after a simple manipulation gives a power of X9^ where x is the dimensionless disorder strength and is a 
small parameter. Since there are no compensating r’s in A^, x can be taken at face value as a small parameter, and 
we expand only to leading order in x- In the dm elements of the diffuson the leading order has a single power of X) 
while all other elements of the D are non-zero at zeroth order in x- Therefore we write 

X,{k + ^ = X,ik) + q-VX,{k) (C5) 

We decompose the diffuson into two parts, a part originating in Xi{k), and another part dg^ originating in the 
gradient VXi{k): 
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This decomposition is physically significant. The s ^ s terms make an exactly null contribution to dg^, so a little 
rearrangement obtains 
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In the last line we have expanded in powers of q, following the same argument as before. In contrast to the previous 
expansion, the first term in the present expansion is matched by a factor of r, so it actually scales as ql without any 
powers of y. Since the denominator subtracts if(s, k) from E{s,k + q) and therefore has no powers of the spin-orbit 
splitting, (ipj can contain no powers of (I -|-4C^)“^. Because of the guarding factor of s multiplying the sum, non-zero 
contributions can be obtained only from a dependence on the spin index s of the Dirac delta function, or of Xi, or of 
VE{s, k). Physically, variations in the Dirac delta function with s reflect a difference between the density of states on 
the spin up s = -1-1 Fermi surface and the density of states on the spin down Fermi surface. Similarly, variations in 
XE{s, k) correspond to differences in the Fermi velocity on the two Fermi surfaces, and variations in Xi to differences 
in the spin-orbit interaction’s angular orientation on the two surfaces. In order to analyze these effects, it is necessary 
to calculate the spin-splitting of the two s = ±1 Fermi surfaces, as a function of the angular variable (pk- It is sufficient 
to obtain them to first order in the ratio of the spin splitting to the Fermi energy, Eprcsn/Ep = xC- After doing so, 
we obtain in the case of the Rashba interaction 
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We have automated this perturbative calculation, and also the calculation of d'”, with mathematica, which facilitates 
the calculation of all matrix elements to high orders. 

We now turn to the other contribution, dg^. Exchanging s,s in the second term and the expanding in powers of 
q-XE (X ql produces 
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The q ■ XXi{k) multiplying the entire expression is proportional to XdP which guarantees a power of y. Therefore all 
other expressions must remain at zeroth order in y, allowing us to neglect the spin splitting effects on the density of 
states, Fermi velocity, and Xi, which are so important for determining dg^. In this approximation the Fermi surface 
becomes a simple circle and the dk turns into a simple angular integral over pk- Moreover, the factor of s multiplying 
everything guarantees that the two s = s contributions sum to zero, leaving us with 

^{-lT?F{k) ■ q)^-^ ^ ^ ( 1 ^ 120 ^ 


For the couplings to in-plane spins, dg^ and dg2, the angular integration selects the M = N term, where N = 1,2,3 
for electrons, polaritons, and holes. For the coupling to out-of-plane spins, dgg, the angular integral selects the M = 2 
term. For iV = 1, 2, 3 the last factor (i=ptT2C)*^ sums to 

(1+4^2)^ ^ (1,2,3(1-4CV3)) (Cll) 

If one performs the remaining angular integral and sums dg^, dg^, then one obtains the charge-spin couplings reported 
in this article. 
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